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Abstract Fast radio bursts show large dispersion measures, much larger than the 
Galactic dispersion measure foreground. Therefore, they evidently have an extragalac¬ 
tic origin. We investigate possible contributions to the dispersion measure from host 
galaxies. We simulate the spatial distribution of fast radio bursts and calculate the dis¬ 
persion measures along the sightlines from fast radio bursts to the edge of host galaxies 
by using the scaled NE2001 model for thermal electron density distributions. We find 
that contributions to the dispersion measure of fast radio bursts from the host galaxy 
follow a skew Gaussian distribution. The peak and the width at half maximum of the 
dispersion measure distribution increase with the inclination angle of a spiral galaxy, 
to large values when the inclination angle is over 70°. The largest dispersion measure 
produced by an edge-on spiral galaxy can reach a few thousand pc cm“^, while the 
dispersion measures from dwarf galaxies and elliptical galaxies have a maximum of 
only a few tens of pc cm“^. Notice, however, that additional dispersion measures of 
tens to hundreds of pc cm“^ can be produced by high density clumps in host galaxies. 
Simulations that include dispersion measure contributions from the Large Magellanic 
Cloud and the Andromeda Galaxy are shown as examples to demonstrate how to ex¬ 
tract the dispersion measure from the intergalactic medium. 
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1 INTRODUCTION 


In the Milky Way, measurements of dispersion measures (DMs) of pulsars have been used to de¬ 
lineate the electron density distribution of the interstellar medium (ISM). Outside our Galaxy, only 
21 extragalactic pulsars in the Magellanic Clouds have DM measurements (iManchester et akl 2006h . 
In recent years, pulsar surveys have discovered a new category of short-duration (a few millisec¬ 
onds) radio bursts which are called fast radio bursts (FRBs). Due to the short duration of FRBs, 
the scale sizes of such exotic astrophysical even ts must be small. Possible origins of FRBs have 
been propose d, such as black hole evapor ation (iKeane et al. 


(lTotanill2013h . binary white-dw arf mergers ( Kashiyama^Lal 


2012h. coalescence of neutron 


stars 


20131) . supermassive rotating neutron 


stars collapsing to black h oles (iFa lcke & Rezzollall2014h . synchrotr on maser emission from rela¬ 
tivistic magnetized shocks (L vubarsk' ^2014t) . magnetar hyperflares (IPopov & Postnovll2010h . and 
near by flaring star s (iLoeb et al. 20141) . Maybe some FRBs are rel ated to gamma-ray bursts (GRBs) 
(e.g. Zhang 2014t). To date, ten such events have been reported ([L orimer et al.ll2007t Keane et al.l 


12012 : Thornton et al.l20f^lSpitler et al.l20T4tlBurke-Spolaor & Bannisteij2014 : |Petroff et al.ll2015 
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iRavi et aDl2015h . Most of these FRBs are found at high Galactic latitudes, and have large DMs up 
to several hundreds pc cm“^, exceeding the possible Galactic foreground DM. Therefore FRBs are 
believed to have an extragalactic origin. 

DMs of a large number of pulsars or FRBs in external galaxies can be used to detect the electron 
density distribu t ion in the diffuse intergalactic med ium (IGM) and ISM in the host galaxies (e.g. 
lDennisonll2014l : Izheng et ^l2014t iHan et al1l2015h . The DM of a given short duration source at 
distance (D, in units of pc) is defined as the integral of the free electron density (rie, in units of 
cm“^) along the line of sight 

DM= / Uedl. (1) 

Jo 

The observed DM, DMobs, of an extragalactic object is the sum of the DM from the host galaxy, 
DMhost, dm in intergalactic space, DMigm, and the foreground DM contributed by our Galaxy, 
DMcai, i-e. 

DMobs = DMhost + DMigm + DMGai- (2) 


The Galactic foreground DM is the integral of electron density along the line of sight from the Sun 
to the furthest reaches of the Milky Way, which can be estimated by different models of free electron 
distribution as discusse d by Schnitze^ ( 201^. The popular thermal electron density model for the 
Milky Way is NE2001 (ICordes & Lazioll2002l) which is mostly built upon pulsar DM measurements. 
However, all known pulsars, except for 21 pulsars in the Large Magellanic Cloud (LMC) and the 
Small Magellanic Cloud, are located inside the Milky Way. Therefore, the electron density distri¬ 
bution beyond pulsars but in the extended Galactic halo cannot be well measured at present, and 
hence the uncertainty of the NE2001 model in the halo region cannot b e assessed. Little is known 
about the DMs contributed by external galaxies. iLorimer et al.l (l2007h estimated a pro bability of 
25% for the DM contribution of 100 pc cm“^ from a host galaxy. Thornton et alJ (12013[) discussed 
the DM contribution of host galaxies and estimated DMhost ~ 100 pc cm“^ in a spiral galaxy 
with a median inclination angle of 60°. The ro ugh distance or redshift estimation of ERB s relies on 
DMigm = DMobs — DMGai — DMhost (see iThornton et al.]l2013l : iLorimer et alJl2007h . with the 
assumption of DMhost = 100 pc cm“^. Although the observed DMs in hig h-redshift galaxies can 
be small due to cosmological time dilation and the associated frequency shift dZhou et al.ll2014l) . the 
local DM contribution from host galaxies to ERBs is very important for calculating the DM from the 
IGM and then estimating their cosmological distances. 

Currently, the origin of ERBs is not known, and the types and morphologies of host galaxies are 
also poorly constrained. In this paper, we apply Monte Carlo simulations to the spatial distribution 
of ERBs in a galaxy, and model the probability distribution of DM from different types of host 
galaxies. The NE2001 model is scaled according to the integrated Ha intensity to represent the 
electron density distributions in host galaxies. We simulate the DM contribution from spiral galaxies 
in Section 2, and discuss DM distributions in other types of galaxies in Section 3. In Section 4 we 
discuss the DM caused by clumps. In Section 5 a specific simulation for the LMC and the Andromeda 
Galaxy (M31) is provided. The conclusions are given in Section 6. 


2 DM CONTRIBUTION FROM SPIRAL GALAXIES 

Most ERBs probably occur in spiral galaxies with locations corresponding to stars or related objects. 
Here we model their spatial distribution and show their probable DM distributions. To specify the 
locations of any modeled ERBs, we use a Cartesian coordinate system (x, y, z) where the center 
of the host galaxy is at the origin. The spatial distribution of ERBs in a galaxy is described by an 
exponential function for height and a Gaussian radial distribution. The surface density of ERBs is 
therefore likely to be large near the galactic center, and decreases with radial distance. The scale 
height of the two-sided exponential function is set to be 330 pc, and the characteristic radius for 
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Fig. 1 DM distributions of FRBs in a spiral galaxy with different inclination angles (0°, 
30°, 60°, 75° and 90°). The solid lines represent the fittings to the distribution with a skew 
Gaussian function. 


the Gaussia n radial density prof ile is set to be 6.5 kpc, similar to that for the simulation of a pulsar 
population dLorimer et al.l 200^ . 

For the Milky Way , the popular model for the density distribution of free electrons is NE2001 
dCordes & Lazioll2002l) . The Galactic DM along a line of sight from the boundary of the Milky Way 
to the Sun can be calculated through the model. For other galaxies the DMs from a host galaxy can 
be calculated along a sightline for the integrals of free electron density from a specific location to the 
near boundary of the galaxy. However, we know little about the electron density distributions in other 
galaxies. For a general spiral galaxy, the size and gas density structure are not very different from 
the Milky Way dBerkhuiisen & FletcheidlTOl^ . We therefore model the electron density distribution 
in a spiral galaxy (size of ~ 30 kpc) using the thin disk, thick disk, spiral arm and Galactic Center 
components of the NE2001 model but ignoring small-scale clumps, voids and the local ISM. 

With this basic model for the electron density distribution, we use a Monte Carlo simulation 
to generate the spatial distribution of ERBs in a spiral galaxy for 10000 locations. When observed 
ERBs come from a distant galaxy outside the Milky Way, all sightlines to these ERBs are toward 
us and nearly parallel. Each spiral galaxy has a specific inclination angle, i, with respect to the 
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Table 1 Fitting parameters of a skew Gaussian function for the DM distributions of FRBs 
in spiral galaxies for ten inclination angles, with the peak DM (column 6) and values for 
the left and right half width at half maximum IFdm-l and Wdm-r- 


i 

(deg) 

(1) 

(pc 

(2) 

LJ 

cm“^) 

(3) 

a. 

(4) 

Hdm-l 

(5) 

DMpeak 
(pc cni~^) 
(6) 

Wdm-r 

(7) 

^^median 
(pc cm“^) 

(8) 

0 

16.7 

38.9 

2.9 

21.6 

35.3 

33.3 

42.5 

10 

16.9 

39.8 

3.0 

21.8 

35.7 

34.2 

43.1 

20 

17.6 

41.9 

3.0 

22.7 

37.3 

36.0 

45.2 

30 

19.3 

45.1 

3.1 

24.3 

40.4 

38.6 

49.0 

40 

23.0 

49.4 

2.7 

28.9 

47.4 

42.2 

55.3 

50 

27.3 

58.5 

2.9 

32.7 

55.4 

50.2 

65.8 

60 

34.4 

75.3 

3.0 

41.6 

70.3 

64.4 

84.0 

70 

49.3 

108.3 

2.9 

60.6 

101.3 

92.8 

121.1 

80 

86.2 

198.8 

2.7 

116.0 

184.1 

170.6 

219.5 

90 

92.7 

527.4 

4.8 

215.3 

293.8 

466.7 

455.2 


sightline. Obviously, FRBs behind the midplane of the disk should show a larger DM than those in 
front of the midplane. We calculate the DM of each FRB by integrating the free electron density in 
the host galaxy along the line of sight from its location to the nearer boundary at 15 kpc from the 
galactic center. Figure [U shows the DM distributions for 10000 assumed FRBs in a spiral galaxy 
with different inclination angles (0°, 30°, 60°, 75° and 90°). The DMs show a peak in the range 
of 30 to 300 pc cm“^, and shift to a large value with increasing inclination angles. The distribution 
follows a skew Gaussian function, defined by 


dN 

dDM 


_ (DM-cF 
= Noe 2..^ 



( 3 ) 


where is the number per DM, and No is the normalization constant; uj and a are the location, 
scale and shape parameters, respectively; t is the variable used in the integral function. The results 
of fitting parameters for the DM distributions are listed in Table [T] The peak DM, DMpeak, and 
the left and right half width at half maximum, VFdm-l and IFdm-r, are given in columns (5) 
- (7) respectively. The peak DMs and the widths at half maximum are plotted in Figure |2] as a 
function of the inclination angle i, both of which increase very quickly at large inclination angles. 
The probability distribution of the inclination angles of spiral galaxies follows the function sin(i), as 
shown in the lower panel in Figured The peak DM exceeds 100 pc cm“^ when the inclination angle 
is over 70°. In a face-on galaxy, when the line of sight passes vertically through the whole disk, the 
DM reaches the maximum of a few tens of pc cm“^, exc ept for the central re gions that have a size of 
^100 pc where the DM can be a few hundred pc cm“^ (iDeneva et al.ll200^ . When the host galaxy 
is edge-on, the sightline to an FRB is at a low height and the very farthest side passes through the 
entire thin disk, so that the DM can be up to 4670 pc cm“^. For the peak DM and the width at half 
maximum, we find empirical formulae to describe DMpeak 33.8 -f 0.94e° °®*, Wdm-l 22.1 -f 
and Wdm-r 39.5 + 0.02e°'^^k We calculate the ensemble distribution marginalized 
over the inclination angle and derive the weighted average of the DM distribution for a spiral galaxy 
to be 142 pc cm“^. If the FRBs follow the distribution of stars with a scale height of young Galactic 
Population I objects 100 pc) (iGilmore & Reidll983h . the DM distribution will be more symmetric 
and have a narrower distribution, but the DM peaks at the same value as shown by the simulations 
above. 
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Fig. 2 The peak DM {squares) and the width at half maximum {the vertical bars) for FRBs 
in spiral galaxies change with the inclination angle i. The dotted line is 100 pc cm“^. The 
lower panel shows the probability distribution of the inclination angles oriented randomly 
in space, which should follow the curve sin(i). 


3 DM CONTRIBUTION FROM OTHER TYPES OE GALAXIES 

FRBs may also happen in dwarf galaxies or elliptical galaxies. The electron density distributions 
of the two types of galaxies are unknown. However, according to their shapes and structures, we 
can model their electron density distributions by scaling the NE2001 model according to the total 
intensity of Ha emission because Ha is a tracer of ionized gas. The emission measure (EM) is the 
pathlength integral of electron density squared, i.e. EM = f nidi. The total Ha luminosity of a 
galaxy is proportional to electron density squared in the volume. 

Based on the galactic structures and their total Ha luminosity, we construct the electron density 
model for dwarf galaxies and elliptical galaxies, and study their DM contribution. Dwarf galaxies 
include irregular galaxies, dwarf spirals and dwarf ellipticals. Here, we are concentrated on young 
and gas-rich Magellanic type galaxies, which are typical dwarf galaxies. The size of Magellanic 
type dwarfs is abo ut 1-5 kpc and the Ha luminosity is 10^^ - 10^® erg s“^ djames et al.ll2004l: 
iKaisin et al.ll2012h . Magellanic type dwarfs are typical irregular galaxies, but they have unconspic- 
uous spiral arms and disks. Eor simplicity, we assume that the size of a dwarf is 3 kpc and the Ha 
luminosity is 10^® erg s“^. The total Ha luminosi ty of the Milky W ay is about 10^® erg s“^, as 
inferred from statistics of Milky Way-like galaxies dJames et al.ll2004h . We scale the Ha luminosity 
of the Milky Way from a diameter of 30 kpc to the size of a dwarf (3 kpc), and set the Ha luminosity 
to 10®^ erg s“^, which is 1/10 of the observed value from a dwarf. The electron density model of a 
dwarf is therefore obtained by scaling the physical size to 1/10 and the amplitude to vTO those of 
the NE2001 model. 

Elliptical galaxies vary greatly in size, from 10 kpc to over 100 kpc. There is very little ISM 
in elliptical galaxies, which results in lo w rates of star format ion. The Ha luminosity of a common 
elliptical galaxy is 10®® - 10^® erg s“^ dKulkarni et al.ll2014ll . The shapes of elliptical galaxies are 
similar to the bulges of spiral galaxies. Here we consider a normal elliptical galaxy with the same size 
as the Milky Way, and its Ha luminosity is 10®® erg s“^. So, in the modeling of electron density in an 
elliptical galaxy, we only use the thick disk and Galactic Center components of the NE2001 model. 
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Fig. 3 The peak value of DM distributions for three types of galaxies as a function of the 
inclination angle i. The width at half maximum is plotted in the vertical bar (scaled to 1/5 
of its original value for clarity). The dotted line is plotted at the value of 100 pc cm“^. The 
dashed line marks the half probability of the inclination angle distribution. 


We then scale the electron density components by reducing the amplitude to i/l/lO according to the 
Ha luminosity. The inclination angle i is dehned by the angle between the line of sight and normal 
direction of the largest plane with the major axis. 

In the Monte Carlo simulations, we scale the spatial distribution of 10 000 FRBs within a galaxy 
accordingly. The DM distributions of FRBs in these galaxies can also be htted with a skew Gaussian 
function as well, but the peak DM and the width at half maximum are much smaller than the corre¬ 
sponding values for spiral galaxies, as shown in Figured FRBs from dwarf galaxies and elliptical 
galaxies have small DMs of a few to a few tens of pc cm“^. The DMs increase slowly from 11- 
12 pc cm“^ at * = 0° to 22 - 24 pc cm“^ at / = 60°, and may reach 100 pc cm“^ when i ^ 90°. 
The ensemble average of the DM distribution for a dwarf galaxy is 45 pc cm“^, and for an elliptical 
galaxy it is 37 pc cm“^. 

4 DM CAUSED BY HIGH DENSITY REGIONS IN A GALAXY 

In the previous sections, only a large-scale diffuse ISM was considered during simulations for DM 
contributions from host galaxies. It is worth noting that there are thousands of small-scale high 
density regions in a galaxy like the Milky Way, where the electron density is strongly enhanced 
relative to the ambient density of the ISM. Though their volume hlling factor is quite small (with the 
maximum of a few tens of pc), the sightline toward an FRB to pass through a clump would produce 
signihcant DM. 

Clumps of ionized gas most likely corres pond to discrete H II regions. To date, more than 2500 
H II regions have been observed in our Galaxy dHou & Hanl2014l) . Their distribution is very inhomo¬ 
geneous and most of them are concentrated in spiral arms. Extragalactic H II regions are also found 
in spi ral galaxies and irregular galaxies, but very few are detected in elliptical galaxies dZhou et al.l 
l2014h . The size of a Galactic H II region ranges from a few percent of a pc (ultracompact) to several 
tens of pc, which is roughly inversely proportional to the electron density in an H II region. The typ¬ 
ical size and average electron density of H II regions detected by the Sino-Germ an 6 cm Polarization 
Survey (e.g. lSun et al.ll20()^ iGao et alJl2O10l : ISun et al.ll201 itlxiao et al.ll201 ih are a few tens of pc 
and a few cm“^, respectively. The DMs caused by the clumps range from tens of pc cm“^ to hun¬ 
dreds of pc cm“^ according to previous studies of pulsars in our Galaxy. Therefore, if an FRB in the 
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Fig. 4 A schematic diagram showing the locations of the LMC and M31 in Galactic coor¬ 
dinates. 


host galaxy happens to take place behind a clump, a large part of the observed DM could come from 
the clump. 


5 APPLICATIONS IN THE LMC AND M3I 

In the near future, a next generation radio telescope, the Square Kil ometre Array (S KA), will reveal 
the electron density distribution in nearby galaxies and the IGM dHan et al.ll2015h . DM measure¬ 
ments of extragalactic pulsars and single pulses are the first powerful probes that can be used to 
detect the extragalactic medium, including ISM in a host galaxy and IGM. A large sample of pulsars 
or FRBs i n nearby galaxies su ch as the LMC and M31 will be detected by this highly-sensitive tele¬ 
scope (see lKeane et alj|2015h . A large sample of Galactic pulsars in the region of sky immediately 
around the direction of the host galaxy can be used to estimate the foreground column density of 
electrons in the Milky Way. After discounting the foreground DM contribution from our Milky Way 
and constraining the local contribution from the host galaxy, the DM from IGM can be derived from 
observations of extragalactic pulses. Here, we apply the simulations for host galaxies of the LMC 

and M31, to demonstrate the detection of intergalactic DM. _ 

The LMC is a satellite of our Milky Way, at a distance of about 50 kpc jPietrzvhski et alJ[2013h 
in the southern hemisphe re of the sky (I = 28 0.5°, b = —32.9°). Its size is about 4.3 kpc, with an 
inclination of ~35° (e.g. lNikolaev et aHl2004t) . At pres ent, 15 pulsars have bee n discovered in the 
LMC, and their DMs have been estimated for 13 pulsars ([Manchester et al.ll2006h. M31 is the larg est 
galaxy in the Local Group, with a distance of approximately 780 kpc ( Stanek & Garnavi^l 199811 in 
the direction (I = 1 21.2°, b = —21.6°). M31 is estimated to have an inclination angle of 77° (e.g. 
IChemin et al.l 2009h . The apparent diameter is about 40 kpc. Their locations relative to the Milky 
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Fig. 5 Simulated DM distributions of extragalactic pulsars (blue) and Galactic pulsars 
(red) within a 10° radius around the LMC on the left and M31 on the right. A small 
sample of pulsars in the LMC and a small sample of foreground Galactic pulsars around 
the cloud are marked with crosses {lower left panel) and presented with as the scaled 
filled area of the histogram {upper left panel). Currently the NE2001 model of electron 
density distribution is not good enough for reconstructing the DMs of foreground pulsar 
data. However, in the future, SKA will discover a much larger sample of Galactic pulsars 
to improve the electron density distribution model and discover many extragalactic pulsars 
(e.g. in M31) that can be used to deduce the DM contribution from host galaxies and the 
IGM. 


Way are illustrated in Figure @1 The simulations mentioned in Sections 2 and 3 are scaled for the 
LMC and M31, respectively, as shown in Figure |5] We simulate foreground Galactic pulsars within 
a 10° radius around the directions of host galaxies. The “observed” DMs from the host galaxies 
in Figure |5] are the sum of DM from the host galaxies, DM in the IGM and the foreground DM 
contributed by our Galaxy in their directions. We add DMigm of 10 pc cm“^ for the LMC, and 150 
pc cm“^ for M31 according to their distance. 

Current data for pulsars in the LMC and in the Milky Way near the direction of LMC are com¬ 
pared with a simulated pulse distribution in the left panel of Figure|5] Obviously, the current NE2001 
model of electron density distribution is not good enough to reconstruct the DMs of farthest Galactic 
pulsars. It is expected that an improved electron density distribution model can be constructed in the 
future. The simulated DMs for pulsars in the LMC are very consistent with the observed DMs of 
LMC pulsars, except for two outliers. 

As pointed out by iHan et all (1201 5h . it is difficult to estimate the dispersion contributed by 
the host galaxy for only an individual extragalactic pulsar or pulse. A large sample of pulsars 
is required for this purpose. As shown in the right panel of Figure |5] the minimum DM of 
these pulsars, min(DMextraPSRs), can be expressed as the upper limit of DMigm plus the fore¬ 
ground Galactic DM, max(DMGaiacPSRs)- The electron density of the IGM can be estimated by 
DMigm = min(DMextraPSRs) - max(DMGaiacPSRs)- When large samples of pulsars or FRBs in 
a host galaxy and also in the Galactic region in the general direction of a host galaxy are discov- 
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ered by SKA in the future, the electron density distribution model for the Milky Way will be greatly 
improved, and the DM contribution from the host galaxy and the IGM can then be deduced. 

6 CONCLUSIONS 

We apply Monte Carlo simulations to generate the spatial distribution of FRBs in a host galaxy. The 
electron density distribution of the host galaxy is modeled by using the scaled model of NE2001. 
The DM distributions of FRBs from the host galaxies are calculated for different types of galaxies, 
which roughly follow a skew Gaussian function. The characteristic parameters of the skew Gaussian 
function increase with inclination angle. FRBs or pulsars in spiral galaxies can have a large DM 
when the inclination angle of the galaxy is over 70° and the sightline of an FRB passes through 
the central region. The largest DM can reach a few thousand pc cm“^ when a galaxy is edge-on. 
The possible DM distributions of pulsars or FRBs in dwarf galaxies and elliptical galaxies are in the 
range of a few to a few tens of pc cm“^. In addition, high density clumps will induce a large DM for 
pulsars or FRBs of tens to hundreds of pc cm“^ if they are behind the clumps. 

In the future when many FRBs or extragalactic pulsars are detected, the statistics of their DM 
distribution can be used to derive the DM contribution by the host galaxy and derive the intergalactic 
DMs after discounting the foreground Galactic DMs. 

Acknowledgements The authors are supported by the National Natural Science Foundation 
of China (Grant Nos. 11473034 and 11503038), and the Strategic Priority Research Program 
“The Emergence of Cosmological Structures” of the Chinese Academy of Sciences (Grant 
No. XDB09010200). 

References 

Berkhuijsen, E. M., & Fletcher, A. 2015, MNRAS, 448, 2469 

Burke-Spolaor, S., & Bannister, K. W. 2014, ApJ, 792, 19 

Chemin, L., Carignan, C., & Foster, T. 2009, ApJ, 705, 1395 

Cordes, J. M., & Lazio, T. J. W. 2002, astro-ph/0207156 

Deneva, J. S., Cordes, J. M., & Lazio, T. J. W. 2009, ApJ, 702, L177 

Dennison, B. 2014, MNRAS, 443, Lll 

Falcke, H., & Rezzolla, L. 2014, A&A, 562, A137 

Gao, X. Y., Reich, W., Han, J. L., et al. 2010, A&A, 515, A64 

Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025 

Han, J., van Straten, W., Lazio, J., et al. 2015, Advancing Astrophysics with the Square Kilometre Array 
(AASKA14), 41 

Hou, L. G., & Han, J. L. 2014, A&A, 569, A125 

James, P. A., Shane, N. S., Beckman, J. E., et al. 2004, A&A, 414, 23 

Kaisin, S. S., Karachentsev, I. D., & Ravindranath, S. 2012, MNRAS, 425, 2083 

Kashiyama, K., loka, K., & Meszaros, P. 2013, ApJ, 776, L39 

Keane, E. R, Stappers, B. W., Kramer, M., & Lyne, A. G. 2012, MNRAS, 425, L71 

Keane, E., Bhattacharyya, B., Kramer, M., et al. 2015, Advancing Astrophysics with the Square Kilometre 
Array (AASKA14), 40 

Kulkami, S., Sahu, D. K., Chaware, L., Chakradhari, N. K., & Pandey, S. K. 2014, New Astron., 30, 51 
Loeb, A., Shvartzvald, Y, & Maoz, D. 2014, MNRAS, 439, L46 

Lorimer, D. R., Bailes, M., McLaughlin, M. A., Narkevic, D. J., & Crawford, F. 2007, Science, 318, 777 
Lorimer, D. R., Faulkner, A. J., Lyne, A. G., et al. 2006, MNRAS, 372, 777 
Lyubarsky, Y. 2014, MNRAS, 442, L9 

Manchester, R. N., Fan, G., Lyne, A. G., Kaspi, V. M., & Crawford, F. 2006, ApJ, 649, 235 


1638 


J. Xu & J. L. Han 


Nikolaev, S., Drake, A. J., Keller, S. C., et al. 2004, ApJ, 601, 260 
Petroff, E., Bailes, M., Barr, E. D., et al. 2015, MNRAS, 447, 246 
Pietrzynski, G., Graczyk, D., Gieren, W., et al. 2013, Nature, 495, 76 

Popov, S. B., & Postnov, K. A. 2010, in Evolution of Cosmic Objects through their Physical Activity, eds. 

H. A. Harutyunian, A. M. Mickaelian, & Y. Terzian, 129 
Ravi, V., Shannon, R. M., & Jameson, A. 2015, ApJ, 799, L5 
Schnitzeler, D. H. E. M. 2012, MNRAS, 427, 664 
Spider, L. G., Cordes, J. M., Hessels, J. W. T., et al. 2014, ApJ, 790, 101 
Stanek, K. Z., & Garnavich, P. M. 1998, ApJ, 503, L131 
Sun, X. H., Han, J. L., Reich, W., et al. 2007, A&A, 463, 993 
Sun, X. H., Reich, W., Han, J. L., et al. 2011, A&A, 527, A74 
Thornton, D., Stappers, B., Bailes, M., et al. 2013, Science, 341, 53 
Totani, T. 2013, PASJ, 65, L12 

Xiao, L., Han, J. L., Reich, W., et al. 2011, A&A, 529, A15 
Zhang, B. 2014, ApJ, 780, L21 

Zheng, Z., Ofek, E. O., Kulkami, S. R., Neill, J. D., & Juric, M. 2014, ApJ, 797, 71 
Zhou, B., Li, X., Wang, T, Ean, Y.-Z., & Wei, D.-M. 2014, Phys. Rev. D, 89, 107303 


